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Abstract 

Modifications to a three-dimensional, implicit, upwind, 
unstructuied-grid Euler code for aeroelastic analysis of com- 
plete aircraft configurations are described. The modifica- 
tions involve the addition of the structural equations of mo- 
tion for their simultaneous time integration with the govern- 
ing flow equations. The paper presents a detailed descrip- 
tion of the time-marching aeroelastic procedure and presents 
comparisons with experimental data to provide an assess- 
ment of the capability. Flutter results are shown for an 
isolated 45° swept-back wing and a supersonic transport 
configuration with a fuselage, clipped delta wing, and two 
identical rearward-mounted nacelles. Comparisons between 
computed and experimental flutter characteristics show good 
agreement, giving confidence in the accuracy of the aeroe- 
lastic capability that was developed. 


Introduction 

In recent years, there has been increased interest in 
the development of aeroelastic analysis methods involving 
computational fluid dynamics techniques. 1 This work in the 
area of computational aeroelasticity has focused on devel- 
oping finite-difference codes for the solution of the tran- 
sonic small-disturbance 2 ’ 3 and full potential equations, 4, 5 
although a growing effort is underway for the solution of 
the Euler and Navier-Stokes equations. 6-17 For example, 
Bendiksen and Kousen 6 presented transonic flutter results 
for two-degree-of-ffeedom (plunging and pitching) airfoils 
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by simultaneously integrating the structural equations of mo- 
tion and the two-dimensional unsteady Euler equations. The 
Euler equations were discretized in space using a finite- 
volume method on a moving mesh and integrated using a 
Runge-Kutta time-stepping scheme. The instantaneous mesh 
was taken to be a superposition of meshes corresponding 
to rigid plunging and pitching of the airfoil. In a follow- 
ing study, Kousen and Bendiksen 7 applied their method 
of Ref. 6 to investigate the nonlinear aeroelastic behavior 
of two-degree-of-ffeedom airfoils at transonic speeds and 
showed that transonic flutter instabilities led to stable limit- 
cycle oscillations involving very large amplitudes. Most 
recently, Bendiksen 8 presented an alternative approach to 
the integration of the structural equations of motion and the 
fluid flow equations. Wu, Kaza, and Sankar 9 integrated in 
time, the unsteady compressible Navier-Stokes equations for 
airfoils undergoing one- and two-degree-of-ffeedom aeroe- 
lastic motions. In Ref. 9, flutter characteristics of airfoils at 
high angles of attack including cases with stall flutter, were 
investigated. The method of Ref. 9 also has been applied 
by Reddy, Srivastava, and Kaza 10 to study the effects of 
rotational flow, viscosity, thickness, and shape on the tran- 
sonic flutter dip phenomena. The study concluded that the 
influence of these effects on flutter, for the cases considered, 
was small near the minimum of the flutter dip, but may be 
large away from the dip. Guruswamy 11-14 demonstrated si- 
multaneous time integration of the three-dimensional Euler 
and Navier-Stokes equations along with the structural equa- 
tions of motion. The inviscid capability first was demon- 
strated in a time-marching flutter analysis performed for a 
rectangular wing with a parabolic-arc airfoil section, and the 
viscous capability was demonstrated later for an aeroelastic 
deformation of a blended wing body to study shock-vortex 
interaction. Finally, Robinson, Batina, and Yang 15 presented 
Euler aeroelastic results for a 45° swept-back wing using a 
deforming mesh capability. A common feature of the above 
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computational aeroelastic methods is that all involved algo- 
rithms that required a structured computational grid. 

As an alternative to structured grid methods, some re- 
cent efforts in computational fluid dynamics are directed to- 
ward solving the governing flow equations using unstruc- 
tured grids. 16-28 These grids typically are constructed from 
triangles in two dimensions, and they consist of an assem- 
blage of tetrahedra in three dimensions. One benefit of 
using these geometric shapes is that they may be oriented 
easily to conform to the geometry being considered, mak- 
ing it possible to represent accurately complicated shapes 
such as multi-element airfoils in two dimensions or com- 
plete aircraft configurations in three dimensions. Applica- 
tion of the unstructured grid methodology has been demon- 
strated for multi-element airfoils 18, 19 and complete aircraft 
configurations. 20, 21 Another benefit of the unstructured grid 
methodology is that the grid data structure simplifies mesh 
refinement in regions of high-flow gradients to resolve the 
physics of the flow more accurately. These adaptive grid 
methods have been demonstrated for steady and unsteady 
flows in two dimensions 22-25 and for complex geometries 
in three dimensions. 26, 27 A disadvantage, however, of the 
unstructured grid methodology is the computational over- 
head of the indirect addressing used to maintain the grid 
data structure. For example, the computational work, for 
steady-state solutions obtained using unstructured grid algo- 
rithms, has been shown to be 2 to 5 times more expensive 
than that required to obtain solutions on a structured grid 
with the same number of cells. 29 

Since there are a number of benefits in using unstruc- 
tured grid methods, it is appropriate to develop computer 
codes for transonic aeroelastic analysis of complete aircraft 
configurations. As a first step, an assessment of the appli- 
cability of the unstructured grid methodology for the aeroe- 
lastic analysis of airfoils was developed and reported by the 
present authors in Ref. 16. Comparisons were made with 
solutions obtained using a structured grid code to determine 
the accuracy of the unstructured grid methodology. The 
unstructured grid capability included a deforming mesh al- 
gorithm to allow the grid to conform to the instantaneous 
position of the moving or deforming airfoil under consider- 
ation. The algorithm is quite general and necessary for the 
treatment of realistic motions encountered during an aeroe- 
lastic calculation. The conclusion in Ref. 16 was that accu- 
rate flutter results could be computed using the unstructured 
grid methodology. In an independent study by Mortchele- 
witcz and Sens, 17 three-dimensional unsteady results were 
presented for a wing undergoing forced harmonic motion. In 
Ref. 17 comparisons of calculated unsteady pressures were 
made with experimental data. The capability included the 
implementation of transpiration boundary conditions that al- 
lowed the mesh to remain fixed for unsteady applications 
where the relative motion of the geometry is assumed small. 

The purpose of the present study is to incorporate 
the aeroelastic analysis procedures of Ref. 16 into a three- 


dimensional, implicit, upwind Euler scheme on unstructured 
deforming meshes and to assess the applicability of the un- 
structured grid methodology for aeroelastic analysis of com- 
plete aircraft configurations. The objectives of the research 
are: (1) to develop a solution algorithm for time accurate 
unsteady flow calculations on a deforming mesh, (2) to im- 
plement the aeroelastic analysis procedures, (3) to compute 
aeroelastic results for an isolated wing and for a complete 
aircraft configuration, and (4) to determine the accuracy of 
the solutions by making comparisons with available exper- 
imental data. The eventual goal is to develop a highly- 
accurate and efficient solution algorithm for the Euler and 
Navier-Stokes equations for aeroelastic analysis of complex 
aircraft configurations. The paper gives a brief description of 
the flow solver used in the current effort, a description of the 
deforming mesh algorithm, and the time-marching aeroelas- 
tic analysis procedures. To demonstrate the time-marching 
aeroelastic procedure that was implemented, flutter results 
are presented for an isolated 45° swept-back wing and a su- 
personic transport configuration with a clipped delta wing 
and two rearward-mounted nacelles. The authors believe 
that these are the first three-dimensional flutter calculations 
obtained using the unstructured-grid methodology. The pa- 
per also presents comparisons between computed and ex- 
perimental flutter characteristics to provide an assessment of 
the accuracy of the capability. 

Upwind-Type Euler Solution Algorithm 

The Euler equations are solved using the three- 
dimensional upwind-type solution algorithm developed by 
Batina. 28 The solution algorithm of Ref. 28 was extended by 
Rausch 30 for time-accurate unsteady flow calculations on a 
deforming mesh and demonstrated for AGARD case 5, pro- 
posed by the AGARD Structures and Materials Panel. 31 The 
algorithm, which is a cell-centered finite-volume scheme, 
uses upwind differencing based on flux-vector splitting, 32 
similar to upwind schemes developed for use on structured 
meshes. The flux-split discretization accounts for the local 
wave-propagation characteristics of the flow and captures 
shock waves sharply with at most one grid point within 
the shock structure. An additional advantage of using flux- 
splitting is that the discretization is naturally dissipative and, 
consequently, does not require additional artificial dissipa- 
tion terms or the adjustment of free parameters to control 
the dissipation. However, in calculations involving a higher- 
order upwind scheme such as this, oscillations in the solu- 
tion near shock waves are expected to occur. To eliminate 
these oscillations, flux limiting usually is required. In the 
present study, a continuously differentiable flux limiter was 
employed. 

The Euler equations are integrated in time using an 
implicit time-integration scheme involving a Gauss-Seidel 
relaxation procedure. 28 The relaxation procedure is imple- 
mented by re-ordering the elements that make up the un- 
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structured mesh from upstream to downstream. The solu- 
tion is obtained by sweeping two times through the mesh 
as dictated by stability considerations. The first sweep is 
performed in the direction from upstream to downstream 
and the second sweep is from downstream to upstream. For 
purely supersonic flows the second sweep is unnecessary. 
This relaxation scheme is stable for large time steps and 
thus allows the selection of the step size based on the tem- 
poral accuracy of the problem being considered, rather than 
on the numerical stability of the algorithm. Consequently, 
very large time steps may be used for rapid convergence to 
steady state, and an appropriate step size may be selected for 
unsteady cases, independent of numerical stability issues. 

Deforming Mesh Algorithm 

A deforming mesh algorithm is used to move the mesh 
for unsteady calculations where the geometry deforms. The 
method, as developed in Ref. 33, models the mesh as a 
spring network where each edge of each tetrahedra repre- 
sents a spring with a stiffness inversely proportional to the 
length of that edge. In this procedure, grid points along the 
outer boundary of the mesh are held fixed, and the instanta- 
neous locations of the points on the wing (inner boundary) 
are specified. For aeroelastic calculations, the position of the 
inner boundary is determined by the structural equations of 
motion. The locations of the interior nodes then are deter- 
mined by solving the static equilibrium equations that result 
from a summation of forces at each node in the x, y, and z 
coordinate directions. The solution of the equilibrium equa- 
tions is approximated by using a predictor-corrector proce- 
dure, which first predicts the new locations of the nodes by 
extrapolation from grids at previous time levels and then 
corrects these locations by using several Jacobi iterations 
of the static equilibrium equations. The predictor-corrector 
procedure is relatively efficient because of the small number 
of Jacobi iterations required to move the mesh. 

Time-Marching Aeroelastic Analysis 

In this section the aeroelastic equations of motion, the 
time-marching solution procedure, and the modal identifica- 
tion technique are described. 

Aeroelastic Equations of Motion 

The aeroelastic equations of motion can be written for 
each mode i as 

+ Ciii + kiqi = Qi ( 1 ) 

where qi is the generalized displacement, m< is the general- 
ized mass, Cj is the generalized damping, k t is the general- 
ized stiffness, and Qi is the generalized force computed by 
integrating the pressure weighted by the mode shapes. These 


equations of motion are derived by assuming that the defor- 
mation of the body under consideration can be described by 
a separation of variables involving the summation of free vi- 
bration modes weighted by generalized displacements. The 
implementation of the equations of motion are slightly differ- 
ent than that of Refs. 3, 15, and 16. In this development the 
generalized aerodynamic forces include modal deflections in 
all three coordinate directions whereas previous implemen- 
tations involved deflections only in the vertical direction. It 
is noted, however, that the results present herein involve 
deflections only in the vertical direction. 

Time-Marching Solution 

The solution procedure for integrating Eq. (1) is similar 
to that described by Edwards et al. 34, 35 A similar formula- 
tion is implemented in the present study for multiple degrees 
of freedom. Here the linear state equations are written as 

ii = Aii + Bui ( 2 ) 

where A and B are coefficient matrices that result from the 
change of variables Xi = [g, q ,} T , and t /,■ is the nondimen- 
sional representation of the generalized force Qi. Equation 
(2) is integrated in time using the modified state-transition 
matrix structural integrator 35 implemented as a predictor- 
corrector procedure. The prediction for i" +1 , x" +1 , is given 
by 

x ? + 1 = fcx? + QB (3u" + <" 1 ) /2 (3) 

where $ is the state-transition matrix, and 0 is the integral 
of the state-transition matrix from time step n to n + 1. 
Then, x" +1 is used to compute the flow field and evaluate 
the nondimensional generalized force u" +1 . These values 
then are used in the corrector step to determine z” +1 , given 
by 

i? +1 = $i? + 0B(^ +1 + u?)/2 (4) 


Modal Identification Technique 

Damping and frequency characteristics of the aeroe- 
lastic responses are estimated from the response curves by 
using the modal identification technique of Bennett and 
Desmarais. 36 The modal estimates are determined by a least 
squares curve fit of the responses of the form 

m 

qi(T) - a 0 + e a ’ T [ajCos(ujT) + bjsin(u/jT )] (5) 

j = i 

i = 1,2,... 

where m is the number of modes. 
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Fig. 1 Planview of 45° swept-back wing model. 



Fig. 2 45° swept-back wing model in the NASA 
Langley Transonic Dynamics Tunnel. 

Results and Discussion 

Flutter results are presented in this section for a 45° 
swept-back wing and for a supersonic transport configura- 
tion with a clipped delta wing and two rearward-mounted 
nacelles. The results are used to assess the time-marching 
aeroelastic capability. The accuracy of the results is deter- 
mined by making comparisons with available experimental 
data. 

45° Swept-back Wing 

To assess the unstructured-grid code for three- 
dimensional aeroelastic applications, a simple well-defined 
wing was selected as a first step towards performing aeroe- 
lastic analyses for complete aircraft configurations. The 
wing that was analyzed was a half-span wind-tunnel-wall- 
mounted model that has a quarter-chord sweep angle of 45°, 
a panel aspect ratio of 1.65, and a taper ratio of 0.66. 37 
The wing is an AGARD standard aeroelastic configuration 



Mode 3, f 3 = 48.35 Hz Mode 4, f 4 = 91.54 Hz 


(a) Oblique projections. 





(b) Deflection contour lines. 


Fig. 3 Natural vibration modes of the 45° swept-back wing. 

which was tested in the Transonic Dynamics Tunnel (TDT) 
at NASA Langley Research Center. A planview of the wing 
is shown in Fig. 1. The wing has a NACA 65A004 air- 
foil section and was constructed of laminated mahogany and 
hence was essentially homogeneous. In order to obtain flut- 
ter for a wide range of flow Mach numbers and densities in 
the TDT, holes were drilled through the wing to reduce its 
stiffness. To maintain the aerodynamic shape of the wing, 
the holes were filled with a rigid foam plastic. A photograph 
of the wing mounted in the TDT is shown in Fig. 2. 

The wing is modeled structurally using the first four 
natural vibration modes which are illustrated in Figs. 3(a) 
and (b). Figure 3(a) shows oblique projections of the natural 
modes, while Fig. 3(b) shows the corresponding deflection 
contours. These modes, which are numbered 1 through 
4, represent first bending, first torsion, second bending, 
and second torsion, respectively, as determined by a finite 
element analysis. The modes have natural frequencies which 
range from 9.6 Hz for the first bending mode to 91.54 Hz 
for the second torsion mode. 

The 45° swept-back wing was modeled using an un- 
structured mesh generated by an advancing front method 
that is part of the VGRID3D 38 software package. The com- 
putational mesh used in the calculations extends two wing 
semispans from the symmetry plane in the span direction. 
Also, the mesh extends ten root chordlengths above/below 
and upstream/downstream of the wing surface to rectangular 
outer boundaries. The meshes on the upper and lower sur- 
faces of the wing are shown in Fig. 4 which indicates that 
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Upper surface Lower surface 


Fig. 4 Surface mesh for the 45° swept-back wing. 



Fig. 5 Partial view of the surface mesh for the symmetry 
plane and the 45° swept-back wing. 

cells have been clustered near the leading edge of the wing. 
The leading edge is a region of large flow gradients, and 
the clustered cells produce a more accurate calculation of 
the leading edge surface pressure. Figure 5 shows a partial 
view of the plane of symmetry and the wing. In this figure 
the mesh along the symmetry plane shows how the cells are 
stretched away from the wing. The complete mesh for the 
wing contains 129,746 tetrahedra and 23,727 nodes. 

Flutter characteristics were calculated for comparison 
with measured values of the flutter speed index and the 
nondimensional flutter frequency at free stream Mach num- 
bers, Moo, of 0.499, 0.678, 0.901, and 0.960 at zero de- 
grees angle of attack. The calculation of each flutter point 
was started by obtaining a steady-rigid solution at the above 
flow conditions. Typically, the next step is to compute a 
static-aeroelastic solution, however, for these steady-state 
flow conditions the wing does not deflect statically since 
the wing is symmetric and at zero degrees angle of attack. 
Therefore, once the steady-rigid solution was obtained for 
each Mach number, a dynamic-aeroelastic calculation was 
started by perturbing the first two structural modes with ini- 
tial velocity conditions. To bracket the flutter point, time- 
marching calculations were performed for several values of 
dynamic pressure, Q, nondimens ionalized by the measured 
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Fig. 6 Comparisons of Euler flutter predictions with 
experimental data for the 45° swept-back wing. 


flutter dynamic pressure, Q exp , including Q/Q exp of 0.7, 
0.8, 0.9, 1.0, and 1.1. 

The aeroelastic responses that result are analyzed using 
the method of Ref. 36 for their damping and frequency com- 
ponents. These components along with their corresponding 
value of dynamic pressure are interpolated to zero damp- 
ing of the dominant flutter mode to obtain the flutter point. 
Figure 6 shows comparisons of computed flutter charac- 
teristics with experimental data. Plots of flutter speed in- 
dex and nondimensional flutter frequency as a function of 
freestream Mach number, are shown in Figs. 6(a) and 6(b), 
respectively. The experimental flutter data defines a typi- 
cal transonic flutter “dip” with the bottom near = 1.0 
for this case. The bottom of the dip in flutter speed index 
(Fig. 6(a)) was defined by the approach to the Moo = 1.072 
flutter point during the wind-tunnel operation. Results from 
the Euler code are presented at the values of Moo at which 
the flutter data was measured. In the Mach number range 
considered in this study (0.499 < Moo ^ 0.96), a conser- 
vative flutter speed was computed at all four Mach num- 
bers in comparison with the experimental data. In general 
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the computed results agree well with experimental data at 
Moo = 0.499 and 0.678 in flutter speed index and in fre- 
quency ratio. Near the transonic flutter dip, however, the 
computed results deviate from the data for the flutter speed 
index but show fair agreement in flutter frequency ratio. The 
flutter results presented above are believed by the authors to 
be the first three-dimensional flutter results obtained using 
the unstructured-grid methodology. 


Supersonic Transport Configuration 

To assess the code for complete aircraft aeroelastic ap- 
plications a calculation was performed for a complex con- 
figuration. This configuration represents an increase in com- 
plexity from that of the 45° swept-back wing from the stand- 
point of the increased complexity of the geometry and the 
natural vibration modes. The configuration analyzed was a 
half-span model of an early supersonic transport (SST) con- 
figuration that was tested in heavy gas in the NASA Langley 
TOT. 39 * 40 This configuration consisted of a rigid fuselage 
and a flexible clipped delta wing with two rearward-mounted 
simulated engine nacelles. A view of the model mounted in 
the TDT is presented in Fig. 7. The wing that is analyzed 
in this paper is that denoted as Wing C in Ref. 39. The wing 
has a leading-edge sweep angle of 50.5° , a panel aspect ratio 



Fig. 7 Supersonic transport model with clipped delta wing 
in the NASA Langley Transonic Dynamics Tunnel. 


of 1.24, and a taper ratio of 0. 142. The airfoil section is a cir- 
cular arc with a maximum thickness-to-chord ratio of 0.03. 
The wing was constructed of a load-carrying aluminum-alloy 
plate structure with cutouts, chemically milled to simulate 
a beam structure and was covered with balsa wood which 
was contoured to the desired airfoil shape. The wing was 
clamped to a relatively rigid mounting block which was at- 



Fig. 8 Deflection contours of the natural vibration modes for the clipped delta wing of the supersonic transport model. 
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tached to a turntable on the tunnel wall. This mounting ar- 
rangement isolated the wing vibrations to the turntable and 
prevented structural interaction between the wing and fuse- 
lage. The model also had two identical slender under-wing 
bodies to simulate engine nacelles. Each nacelle consisted 
of a cylindrical centerbody with an ogive nose section and 
a conical tail fairing. The total mass of the nacelles was 
about the same as the total mass of the wing. The fuselage 
fairing was a half body of revolution that was extended from 
the tunnel-wall to ensure that the wing root was outside the 
tunnel-wall boundary layer (Fig. 7). 

Nine natural vibration modes and their associated gen- 
eralized masses were measured. 39 Deflection contours of 
these wing modes are shown in Fig. 8. These modes have 
natural frequencies that range from 7.8 Hz for mode 1 to 
58.1 Hz for mode 9. The nacelle masses have a large effect 
on the mode shapes, as shown in Fig. 8, particularly in the 
inboard region of the wing. 

The SST configuration also was modeled using the 
VGRID3D mesh generation package. The computational 
mesh extends two wing semispans from the symmetry plane 
in the span direction. Also, the mesh extends ten root 
chordlengths above/below and upstream/downstream of the 
wing surface to rectangular outer boundaries. The top, bot- 
tom, and side views of the surface mesh of the configura- 
tion are shown in Fig. 9. The upper and lower views of the 
surface mesh show that cells have been clustered near the 
wing tip and around the nacelles. The side view shows that 


the wing is placed below the centerline of the fuselage. The 
complete mesh for the supersonic transport contains 323,818 
tetrahedra and 59,429 nodes. 

The measured natural vibration mode shapes are inter- 
polated to the surface mesh of the configuration of Fig. 9. 
The interpolated mode shapes of the configuration are shown 
in Fig. 10 along with the corresponding natural frequen- 
cies. Figure 10 illustrates the relative vertical motion of 
the clipped delta wing and nacelles with respect to the rigid 
fuselage. Similar to the model in the wind-tunnel test, the 
fuselage was rigid as shown in Fig. 10. It should be noted, 
however, that the time-marching aeroelastic capability al- 
lows for general motion of the complete configuration and 
is not restricted to simple wing deflections. 

A calculation was performed for the SST configuration 
at Moo = 0.907 and zero degrees angle of attack. The 
aeroelastic calculation was performed by first obtaining a 
steady-rigid solution at these flow conditions. Next a static- 
aeroelastic solution, during which the wing was allowed 
to deform due to the aerodynamic loads caused by the 
nonsymmetric geometry, was computed. To allow rapid 
convergence of the wing to its static deformed shape and 
to prevent the wing from oscillating, structural damping 
was added. Finally a dynamic-aeroelastic calculation was 
started from the static-aeroelastic solution by perturbing the 
first three structural modes with initial velocity conditions. 
Figure 1 1 shows the resulting generalized displacements for 
the first three structural modes, where a value of dynamic 




Fig. 9 Surface mesh of the supersonic transport configuration with a clipped delta wing and two rearward-mounted nacelles. 
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Mode 4, f 4 = 25.4 Hz Mode 5, f 5 = 38.2 Hz Mode 6, f 6 = 43.3 Hz 



Mode 7, f 7 = 45.9 Hz Mode 8, f 8 = 48.2 Hz Mode 9, f 9 = 58.1 Hz 

Fig. 10 Natural vibration mode shapes of the supersonic transport configuration with a clipped delta wing and two 
rearward-mounted nacelles. 

pressure that was found experimentally to correspond to transient is near the flutter point Therefore, the computed 

flutter was used. The first three component modes of value of dynamic pressure is in good agreement with the 

the second generalized displacement are shown in Fig. 12. value found experimentally to correspond to flutter. The 

Damping and frequency estimates of the three modes are calculated flutter frequency was found to be 10 Hz which 

listed in Table 1. The near neutrally stable Mode 1 of the compares with the experimental flutter frequency of 11 Hz. 

generalized displacement indicates the computed aeroelastic 
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Fig. 11 Generalized displacements at Q/Q exp = 1.0 
for the supersonic transport configuration at 
Moo = 0.907 and zero degrees angle of attack. 


Concluding Remarks 

Modifications to a three-dimensional, implicit, upwind 
Euler code based on unstructured grids for the aeroclastic 
analysis of complete aircraft configurations were described. 
The modifications involved the addition of the structural 
equations of motion for their simultaneous time integration 
with the governing flow equations. The flow solver of the 
Euler code, which is a cell-centered finite-volume scheme, 
uses upwind differencing based on flux-vector splitting and 
involves an implicit time-integration scheme which uses a 
Gauss-Seidel relaxation procedure. The code also includes 
a deforming mesh algorithm that is capable of moving the 
mesh for general aeroelastic motions of complete aircraft 
configurations. Flutter results w ere presented for an isolated 
45° swept-back wing and a supersonic transport configura- 
tion with a clipped delta wing and two rearward-mounted na- 
celles to assess the time-marching aeroelastic procedure that 
was implemented. Comparisons show good agreement be- 
tween computed and experimental flutter characteristics, giv- 
ing confidence in the accuracy of the aeroelastic capability. 
The authors believe these Me the first three-dimensional flut- 
ter calculations obtained using the unstructured-grid method- 
ology. 


Table 1 Component mode damping and frequencies 
of the aeroelastic system at Q/Q exp - 1.0 
for the supersonic transport configuration at 
Moo = 0.907 and zero degrees angle of attack. 


Mode 1 

Mode 2 

Mode 3 

c 

u (Hz) 

C 

v (Hz) 

c 

u (Hz) 

0.002 

10.06 

0.082 

16.04 

0.027 

24.08 



Fig. 12 First three component modes of the second 
generalized displacement at Q/Q exp =1.0 
for the supersonic transport configuration at 
Moo = 0.907 and zero degrees angle of attack. 
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